Investigation of Single Boron Acceptors at the Cleaved Si:B (111) Surface 
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The cleaved and (2 x 1) reconstructed (111) surface of p-type Si is investigated by scanning tunneling 
microscopy (STM). Single B acceptors are identified due to their characteristic voltage-dependent 
contrast which is explained by a local energetic shift of the electronic density of states caused by the 
Coulomb potential of the negatively charged acceptor. In addition, detailed analysis of the STM 
images shows that apparently one orbital is missing at the B site at sample voltages of 0.4 — 0.6 V, 
corresponding to the absence of a localized dangling-bond state. Scanning tunneling spectroscopy 
(STS) confirms a strongly altered density of states at the B atom due to the different electronic 
structure of B compared to Si. 



PACS numbers: 6I.16.Ch, 68.35. Bs, 73.20.-r 
I. INTRODUCTION 

In recent years much interest has been devoted to iden- 
tify individual dopant atoms at the surface of doped semi- 
conductors by scanning tunnehng microscopy (STM) . For 
instance, on GaAs (110) substitutional Si donors (Sica) 
and Be or Zn acceptors in the top few surface layers ap- 
pear as protrusions due to the local change of the tip- 
induced band-bending arising from the Coulomb poten- 
tial of the ionized dopant |I|g] . At low temperatures and 
negative sample bias circularly modulated structures in 
the topographic image have been interpreted as Friedel 
oscillations around the ionized donors induced by the ac- 
cumulated electrons in the space-charge region H. Fur- 
thermore, in n-dopcd InAs the scattering states of ionized 
dopants at low temperatures have been explored only re- 
cently |. 

In addition to the delocalized features in the subsur- 
face region, localized features have also been found which 
were attributed to a local change of the electronic struc- 
ture around dopants at the GaAs (110) surface ||l[. Re- 
markably, ah initio calculations predict that the addi- 
tional electron of the substitutional Sioa atom is trapped 
by a localized midgap level due to a local modification 
of the electronic structure around the Si atom g] . Fur- 
thermore, different kinds of dopant-induced features and 
surface defects have been distinguished by voltage depen- 
dent imaging of the occupied and unoccupied electronic 
states g. 

Apart from studies of the segregation or adsorption 
of group III or group IV elements on the Si (111) sur- 
face, most of the work on the local electronic structure 
around individual dopant atoms focussed on doped III-V 
semiconductors rather than on elemental semiconductors 
like Si, although these surfaces have been investigated in 
great detail by STM and scanning tunneling spectroscopy 
(STS) 0. The question of how the dopant atoms are 
spatially distributed is of crucial interest for the investi- 
gation of the metal-insulator transition in Si doped with 
phosphorus and/or boron pi. In a previous report we 



have shown that individual P donors at the (2x1) re- 
constructed (111) surface of n-doped Si can be identified 
by voltage dependent imaging and STS at room tempera- 
ture P,EO[ . Electronic surface states energetically located 
in the bulk band gap pin the Fermi level at 0.4 eV above 
the top of the valence band at the surface. Consequently, 
the bands are bent upward towards the surface in con- 
trast to GaAs (110) where the bands remain flat at zero 
bias. In Si:P the Coulomb potential around the ionized 
donor causes a local down-shift of the electronic density 
of states (DOS) and the donor appears as a protrusion 
at positive sample bias and as an indentation at negative 
bias. 

In order to further investigate the change of the local 
electronic structure around dopant atoms at the Si (111) 
surface we performed STM and STS measurements on 
cleaved p-doped Si (Si:B) at room temperature. Previous 
studies on the adsorption of B on Si (111) revealed dif- 
ferent stages of B incorporation in the surface depending 
on coverage and thermal treatment ||ll|,|l^ in contrast to 
other group III adsorbates p3|jl^ ]. For the cleaved Si:B 
(111) surface one might naively expect that apart from 
the change of the sign of the Coulomb potential the influ- 
ence of the negatively ionized B acceptor on the electronic 
structure can be explained in a way similar to Si:P. How- 
ever, we will show that unlike in Si:P, the DOS in Si:B 
is strongly modified at the acceptor site possibly due to 
the different electronic configuration of substitutional B 
in Si compared to P. 



II. EXPERIMENTAL 

Measurements were performed with an Omicron STM 
in ultra-high vacuum (UHV) at room temperature. STM 
tips were prepared from electro-chemically etched tung- 
sten wire and further cleaned in UHV by repeated cycles 
of annealing and consecutive Ar+ sputtering. Samples 
(0.3 X 4 



X 10 mm ) with boron concentrations Na 

^19 



7-10 cm •^ and 4.5-10 cm ■^ were cut from Czochralski- 
grown single crystals [la and cleaved in situ to expose 



the (111) surface to the tip. STM images were acquired 
directly after cleavage without further heat treatment of 
the samples to maintain the original dopant distribution 
at the surface. Images were taken in the constant-current 
mode with the voltage applied to the sample and the tip 
grounded. Hence, at positive voltages unoccupied elec- 
tronic states of the sample are imaged whereas at nega- 
tive voltages occupied states are imaged, implying that 
the DOS of the tip varies smoothly with energy. Scans of 
opposite polarity were acquired quasi-simultaneously by 
scanning each line forward and backward with reversed 
polarities. 



III. RESULTS AND DISCUSSION 

A. Identification and Distribution of B Acceptors 

Fig. H shows STM images of the cleaved Si(lll) sur- 
face at -1-1.2 V and —1.2 V respectively. The bright rows 
are characteristic for the (2 x 1) reconstruction of the 
cleaved (111) surface as investigated in detail by Feen- 
stra et al. ||l^. The reconstruction is explained by a 
revised 7r-bonded chain model ||l^j|l3l where the Pz or- 
bitals representing 'dangling bonds' (DB) retain their lo- 
cal character and are therefore ideally suited for imaging 
with an STM. Fig. S shows a schematic sketch of this 
reconstruction where atoms marked 3 and 4 are fourfold 
coordinated and lowered, while the outer atoms 1 and 
2 have one DB each which is mainly occupied at atom 
1 and mainly unoccupied at atom 2 pSJ. In the STM 
only the elevated atoms 1 and 2 are imaged. At neg- 
ative sample voltage electrons tunnel out of the mainly 
occupied orbitals at position 1 which therefore appear as 
bright chains in the image, separated by dark stripes cor- 
responding to atoms at positions 3 and 4. At low positive 
voltages the electrons are able to tunnel into the mainly 
unoccupied states located at atom 2 which therefore ap- 
pear as bright chains. The small shift of the rows along 
the [211] direction is verified by comparing images of op- 
posite polarity. At even higher positive voltages elec- 
trons tunnel into unoccupied states that are located at 
the bonds between atoms 1 and 2 so that the rows appear 
as a zigzag chain. Several types of defects are observed. 
Those marked by arrows have a characteristic voltage- 
dependent contrast. These defects appear as protrusions 
(bright) at —1.2 V in contrast to images taken at -fl.2 
V where they appear as indentations (dark). Hence, the 
contrast is exactly inverted with respect to Si:P. These 
characteristic defects are observed for both investigated 
boron concentrations. We count a total of 67 defects in a 
surface area of 21323 nm^ for Na = 7- lO^^cm'^ and 106 
defects in 8357 nm^ for Na = 4.5 • lO^^cm-^. The values 
correspond to dopant densities of 3.14 • 10^'^nm^^ and 
1.26 • 10~^nm~^, in good agreement with the respective 
surface densities of 2.2 • 10~'^nm~^ and 1.4 • 10~^nm~^ 




FIG. 1. STM images (30 nm x 30 ran) of the cleaved Si:B 
(111) surface {Na = 7 • lO^^cm"^) at negative and positive 
sample voltage U , tunneling current / = 0.7 uA. White arrows 
indicate B-induced features. 
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FIG. 2. Schematic sketch of the 2x1 reconstructed (111) 
surface. The lobes indicate the occupied or unoccupied dan- 
gling- bond states at atom 1 or 2, respectively. 

obtained from the bulk concentrations, where we as- 
sume that only atoms in the outermost layer (atoms 1 
to 4, Fig. H) give rise to a contrast in the image. We 
therefore ascribe these features as being due to single B 
acceptors in the surface layer. It is reassuring that struc- 
tures with this specific voltage-dependent contrast have 
not been found on the previously investigated Si:P (111) 
surface Hjl^. Vice versa, the characteristic voltage- 
dependent contrast ascribed to P donors on Si:P has not 
been found on the Si:B surface investigated here. 

The identification of these defects as individual B 
dopants allows to check whether the dopant distribution 
in Si:B is random. The probability of finding the nearest 
neighbor B atom at a distance r from a given atom at 
the surface is /(r) = 27rrpe~''^'' , assuming a Poisson dis- 



tribution of dopants with the dopant surface-density p. 
For the determination of the nearest neighbor distances, 
B atoms being nearer to the image border than to any 
other B atom have been discarded. The distances r have 
been grouped into 0.665 nm wide intervals corresponding 
to the unit cell dimension along [211]. Fig. shows a his- 
togram together with /(r) normalized to the area under 
the histogram. The behavior does not change consider- 
ably when the intervals are shifted by the row-to-row dis- 
tance of 0.332 nm. We find reasonably good agreement 
between the statistical and experimental distribution, al- 
though there seems to be a cut-off at low r similar to 
what has been observed for Si:P. The shortest distance 
between two B atoms was found to be 1.67 nm. This is 
in reasonable agreement with the distance of 2.3 nm that 
follows from the solubihty limit 1.2 at% (T < 1400 K) 
of B in Si for a random distribution [|l^. The reduced 
counts at longer distances are likely to be due to the re- 
stricted size of the STM images. For this reason, the 
statistics has been carried out for the higher concentra- 
tion of Na = 4.5 • lO^^cm^'^ only, because most of the im- 
ages taken on the sample with Na = 7- lO^^cm"'^ showed 
only one or two B defects. We conclude that B acceptors 
in Si:B are distributed randomly with a short-distance 
cut-off similar to P donors in Si:P. 
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FIG. 3. Distribution of B nearest neighbors as determined 
from STM images (histogram) and /(r) derived for a random 
arrangement of B atoms (solid line) . 



B. Electronic Structure of Single B Acceptors 

We now discuss the voltage-dependent contrast in more 
detail. Fig. ^ shows STM images of a B-induced feature 
for different voltages. At all negative voltages boron ap- 
pears as an isotropic shallow protrusion with a diameter 
of ~ 2 nm. In contrast, at positive voltages B appears as 
a needle-shaped protrusion at -1-0.4 V but can hardly be 
distinguished from Si at -1-0.9 V. Further increase of the 
voltage leads to a reversal of the image contrast compared 
to -1-0.4 V and the acceptor appears as an indentation. 
This change in contrast with voltage is explained by a 
local change of the surface DOS around the acceptor. 




FIG. 4. STM images (7x7 nm^) of the Si:B (111) surface 
{Na = 7- 10^**cm"^) at different sample voltages, / = 0.7 nA. 

Fig. H shows the surface DOS for the p-type cleaved Si 
(111) surface together with the near-surface bulk energy 
bands as inferred from photoemission and STS experi- 
ments pmEOp^]. In p-type Si, the acceptor level at 45 
meV above the valence band edge is occupied at room 
temperature and the boron atom is negatively charged. 
This occurs for B acceptors at the surface as well, due to 
the low hole binding energy of 25 meV estimated for high 
carrier concentrations P2| . In addition, due to the pres- 
ence of electronic surface states in the bulk band gap (Fig. 
|5|, solid line), holes (majority carriers) are accumulated 
at the surface and the positive charge is compensated by 
the formation of a hole depletion layer of depth d. The 
negative space charge of the depletion layer gives rise to 
a downward band bending towards the surface. Further- 
more, the Fermi level is pinned 0.4 eV above the valence 
band edge at the surface almost independent of dopant 
concentration for moderately doped Si |20 2^. For a 



boron concentration of Na = 4.5 • 10 cm"'' and a band 
bending of —0.4 V a depletion-layer depth d = 3.4 nm 
and a surface charge corresponding to 1.5 x 10^'^ ecm~^ 
are estimated [g3[. This charge corresponds to a large 
hole density of 2.6 x 10^° cm~'^ near the surface owing 
to the fact that the surface states decay into the volume 
with a short decay length of 1 - 2 lattice constants [p4|. 
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FIG. 5. Surface density of states and energy bands for 
p-doped Si (111) at the Si position (solid line) and at the 
acceptor position (dashed line). Hatched areas indicate occu- 
pied electronic states of Si. 

Adopting the simple description which worked success- 
fully for Si:P, the surface DOS at the location of the boron 
atom is energetically shifted upward due to the screened 
negative Coulomb potential of the B acceptor (Fig. |^, 
dashed line). For the occupied states below Ep this up- 
ward shift leads to a local increase of the surface DOS 
near the Fermi level. At negative voltages, electrons near 
the Fermi level of the sample have the highest transmis- 
sion factor for tunneling into unoccupied states of the tip. 
Therefore, at the site of the B atom more electrons are 
able to tunnel into the tip compared to Si sites. Thus, B 
appears as a protrusion for all negative voltages. 

At positive voltages, electrons at energies near the 
Fermi level of the tip have the highest probability for 
tunneling into unoccupied surface states of the sample. 
At -1-0.4 V there are more states available at the position 
of B which therefore appears as a bright protrusion. At 
-1-0.9 V the surface DOS of B and Si are alike so that B 
and Si can hardly be distinguished. At even higher volt- 
ages (-1-1.2 V) there are more surface states on Si, lead- 
ing to a higher tunneling probability at the Si position 
than at the B position and B appears as an indentation. 
Hence, in this simple picture of a locally shifted DOS 
the contrast changes from bright to dark with increasing 
positive voltage as experimentally observed. 

A more detailed analysis of the image taken at -1-0.4 V 
shows that apparently one unoccupied orbital in the row 
containing the defect is missing, compared to the image 
for U = —0.4 V (Fig. g). The same behavior was found 
for ±0.6 V bias. The fact that the orbital is present for 
higher voltages confirms that the feature is not due to a 
surface vacancy. The missing orbital is better seen in the 
cross-section lines where each maximum corresponds to 
an unoccupied orbital (Fig. ph. At the position of the 
B acceptor a minimum appears at a position where ac- 
tually a maximum is expected. Counting the number of 
maxima along the two lines B, B', i.e. along the [Oil] di- 



rection, one maximum is missing in the line taken across 
the defect (B) compared to a scan away from the acceptor 
(B'). We stress that this behavior was found for about 50 
% of all B-induced features while 50 % appeared regular. 
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FIG. 6. STM images (7 nm x 7 nm) of the boron-induced 
feature for U = ±0.4 V (cf. Fig. q) and height profiles Az{r) 
along the indicated lines. In the scan along B one orbital is 
missing at the position of boron compared to the scan along 
B'. 

In Si:P the appearance of an additional orbital at the 
position of the P donor at U = —0.4 V ||9|,pi| was ex- 
plained as being due to Si surface states just above Ep, 
which arc usually unoccupied but are shifted to below 
the Fermi level by the Coulomb potential of the positively 
charged P donor. These states refer to the dangling-bond 
orbital at atom 2 (cf. Fig. 0). Therefore, this orbital be- 
comes occupied and appears as an additional "atom" at 
negative voltages. As is apparent from Fig. 5, a cor- 
responding explanation invoking simply a local shift of 
the surface DOS cannot be found to account for a miss- 
ing orbital at the acceptor site in Si:B. We believe that 
this is due to different substitutional positions of boron 
at the surface. Rather, a qualitative change of the sur- 
face DOS at the B site must be invoked. For instance, 
the trivalent B located at position 1 or 2 would partici- 
pate in three covalent bonds although the preferred flat 
sp^-hybrid configuration cannot be realized. In contrast, 
B at position 3 or 4 has to satisfy four bonds making a 
charge transfer from adjacent Si atoms very likely. Such 
an effect of the B substitutional site on the electronic 
structure has been reported for B deposited on Si (111) 
where B does indeed occupy immediate subsurface posi- 
tions |l^Jl^. In addition, the Si-B bond length is about 
12 % shorter than the Si-Si bond length mainly due to 
the smaller covalent radius of B compared to Si iQ . This 
will lead to a relaxation of the surface structure around 
B which could be connected with a charge transfer be- 



tween Si and B and/or a reformation of the unoccupied 
DB at/nearest to the B site towards a more back-bonded 
type. A missing orbital in the STM image at +0.4 V 
could either be explained by a complete occupation of 
the orbital at position 2, which then would not be im- 
aged for U > 0, or by a reformation of the DB towards a 
more back-bonded character. 
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FIG. 7. Normalized tunneling spectra taken away from and 
on top of the B-induced feature, respectively. Different sym- 
bols indicate different tunneling spectra calculated from an 
averaged set of several I{U) curves. 

Because only atoms 1 or 2 are imaged the effect of B 
at position 3 or 4 on the image contrast is presumably 
weaker because the acceptor at position 3 or 4 can only 
be "imaged" via an electronic interaction to Si atoms at 
positions 1 or 2. Hence, we conclude that the effect of 
a missing orbital is presumably due to B located at po- 
sitions 1 or 2. This is in good agreement with the fact 
that we found a missing orbital at roughly 50 % of the 
B sites. We emphasize that the topological structure of 
the B-induced feature clearly requires theoretical calcu- 
lations. Moreover, a more sophisticated model has to in- 
clude many-body effects which are important for highly 
localized electronic states WM . 

The local change of the DOS at the acceptor site is con- 
firmed by scanning tunneling spectroscopy (STS). During 
I{U) data acquisition the tip was retracted from the sur- 
face by 0.05 nm/V to compensate for the exponential in- 
crease of current with voltage. Several I(U) curves were 
averaged to reduce the scatter of the data. As demon- 
strated earlier, the logarithmic derivative {dI/dU)/{I/U) 
is independent of the distance and is related to the sur- 
face DOS with U = corresponding to Ep pM ■ Due to 



the numerical evaluation of {dI/dU)/{I/U) data points 
around U = were omitted because division by zero leads 
to unreasonable results, while {dl / dU) / {I /U)\u=o — 1. 
Fig. M shows spectra taken above a Si atom away from the 
acceptor which exhibits peaks at U = —0.9, -1-0.5, -1-1.4 V 
in agreement with previous results for p-doped Si where 
similar peaks were observed at somewhat lower voltages 
U = -1.0, -0.35, -h0.17, -Hl.25 V due to the lower doping 



level 1 16 compared to the present samples. These peaks 
can be attributed to the surface DOS of Si observed in 
photoemission experiments [ p0|j2l] . A peak expected at 
U = —0.35 V corresponding to the occupied DB states 
could not be resolved due to the numerical difficulties 
around U = 0. 

In comparison, the spectrum taken above the acceptor 
is strongly modified, apart from the peak at C/ = -I-1.45V 
which is slightly shifted by AC/ < 0.1 V compared to 
the peak on Si. The most striking result is the strong 
reduction of the large peak observed on Si at -1-0.5 V. The 
corresponding unoccupied electronic states are located at 
the DB orbitals at position 2. Strong modifications of the 
STS spectra have also been reported for B deposited on 
Si (111) |l|,|lj]. In the present case the reduction of the 
surface DOS above Ep confirms the absence of a localized 
DB state at the B site and the effect of a missing orbital 
in the STM image at [/ = 0.4 — 0.6 V. Furthermore, the 
peak at —0.9 V representing the bottom of the occupied 
surface band is considerably shifted upward by « 0.5 V 
at the B site. Since at negative bias, states near Ep of 
the sample dominate the tunneling current this conffrms 
the bright image contrast at the B site observed for all 
negative voltages. 

A further question concerns the apparent anisotropic 
shape of the B defect. As mentioned above, at —0.4 V the 
extension of the B induced feature seems to be isotropic 
while it appears strongly elongated along the [Oil] direc- 
tion at positive voltage. Such an effect could be due to 
a different electronic interaction between Si and B along 
the vr-bonded chains compared to the perpendicular di- 
rection and needs further investigation. We mention that 
the P-induced feature on Si:P appeared to be isotropic, 
although its overall extension was smaller. An explana- 
tion for the larger extent of the B-derived features could 
be a different distance between the tip and the surface 
compared to Si:P caused by a higher tunneling current 
of 0.7 nA compared to 0.3 nA in the latter case. 



IV. SUMMARY 

The Si (111) (2 x 1) surface of cleaved B-doped single 
crystals has been investigated by STM. The individual 
B acceptors have been identified by a voltage-dependent 
contrast and have been further characterized by STS. 
Similarly to Si:P the general change in contrast is gov- 
erned by the local shift of the surface DOS due to the 



Coulomb potential of the charged dopant atom. In addi- 
tion, the local change of the atomic and electronic struc- 
ture around the B acceptor gives rise to a missing orbital 
observed in 50 % of the images taken at U = -1-0.4 V. 
These acceptors are presumably located at the upper po- 
sitions 1 and 2 of the (2x1) reconstructed surface. The 
unoccupied DOS of B at the surface as determined by 
STS is strongly reduced around 0.4 eV, in contrast to 
Si:P where such drastic changes are not observed. This 
shows that the electronic configuration of substitutional 
dopants has a decisive influence on the local electronic 
structure. 

The fact that in Si:B and Si:P the surface-densities 
of dopants estimated from the STM data are in good 
agreement with the surface densities derived from the 
bulk concentrations strongly confirms the assumption 
that only dopants in the outermost layer are imaged. 
This is in strong contrast to GaAs, where features aris- 
ing from dopants in several subsurface layers are observed 
[0,||. Presumably, the screening of dopants is very dif- 
ferent in doped Si and needs to be investigated further. 
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